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Abstract 

Astrophysical sources of ultrahigh energy neutrinos yield tau neu- 
trino fluxes due to neutrino oscillations. We study in detail the con- 
tribution of tau neutrinos with energies above 10^ GeV relative to the 
contribution of the other flavors. We consider several different ini- 
tial neutrino fluxes and include tau neutrino regeneration in transit 
through the Earth and energy loss of charged leptons. We discuss sig- 
nals of tau neutrinos in detectors such as IceCube, RICE and ANITA. 



1 Introduction 



Very high energy neutrinos can be unique probes of both astrophysics and 
particle physics. They point back to their sources, escape from the most 
extreme environments and have energies much higher than those available in 
colliders. 

The SuperKamiokande experimental data on atmospheric neutrinos 
shows evidence of nearly maximal ^ Vr mixing. As a result, astrophys- 
ical sources of z/^ become sources of and Vr in equal proportions after 
oscillations over astronomical distances |2|. We do not differentiate between 
neutrinos and antineutrinos, since they can not be distinguished experimen- 
tally. The neutral current cross sections for and are identical, and 
above ~ 1 TeV the charged current cross sections are identical as well. Nev- 
ertheless, signals from Ur have the potential to contribute differently from 
z/^'s because their propagation through matter is different. Tau neutrinos 
produce tau leptons via charged current interactions in the Earth. Having a 
short lifetime, the tau leptons decay producing Vj. which then interact and 
produce r, resulting in a cascade that produces z/^ and r with energies lower 
than the original flux OE]. The leptonic tau decays also produce secondary 

and Ve neutrinos All neutrinos from this cascade can then interact 
in the detector. The decays of taus in the detector also contribute to the 
signal E] ■ For muons, the electromagnetic energy loss coupled with the 
long muon lifetime make regeneration of from muon decay negligible and 
high energy z/^'s get large attenuation as they propagate through the Earth. 

Signals of neutrino interactions in the rock below the ice or in ice depend 
on the energy and flavor of the neutrino. Muon neutrino charged current 
(CC) conversions to muons are noted by the Cherenkov signal of upward 
going muons in a detector such as IceCube p|. High energy electromagnetic 
showers from z/g — > e CC interactions produce Cherenkov radiation which 
is coherent for radio wavelengths [H]. The Radio Ice Cherenkov Experiment 
(RICE) has put limits on incident isotropic electron neutrino fluxes which 
produce downward-going electromagnetic showers ^j. The Antarctic Impul- 
sive Transient Antenna (ANITA) also uses the ice as a neutrino converter 
|lUj . These balloon missions will monitor the ice sheet for refracted radio 
frequency signals with an effective telescope area of IM km^. All flavors of 
neutrinos produce hadronic showers. In addition, tau decays contribute to 
both electromagnetic and hadronic showers that could be detected by Ice- 
Cube, RICE or ANITA. 
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In this paper, we investigate the effect of z/^ regeneration from tau de- 
cays and tau energy loss for neutrinos with energies above 10^ GeV, with a 
particular interest in the higher energy range relevant for RICE and ANITA. 
Attenuation shadows most of the upward-going solid angle for neutrinos, so 
we concentrate on incident neutrinos which are nearly horizontal or slightly 
upward-going. Considering the correct propagation of tau neutrinos is im- 
portant for a number of reasons. The effective volume for the detection of 
is larger than for the other flavors due to the regeneration effects. The initial 
interaction happens far outside the detector, but the neutrinos (taus) pro- 
duced in the neutrino interaction-tau decay cascade interact (decay) inside 
the detector. This can lead to enhanced event rates for particular energies 
and trajectories. The regenerated neutrinos and taus contribute, however, 
at lower energies than the initial ones, so a detailed quantitative analysis is 
necessary to understand where the regeneration effects are significant. Be- 
sides the possible enhancement in the rates, a good discrimination between 
all the different neutrino flavors is important. The flavor composition of the 
detected neutrino fluxes could provide a better understanding of the sources 
of such neutrinos and neutrino properties. The separate identiflcation of each 
flavor could be important for this purpose [TT]. 

We illustrate our results with a variety of fluxes. We concentrate our 
analysis on the flux of GZK neutrinos ^21- These are produced in the de- 
cay of pions from the interaction of cosmic ray protons with the background 
microwave photons, so they are a "guaranteed" source. In addition, they 
are expected to dominate the neutrino flux at energies for which the radio 
Cherenkov detection methods are most sensitive if the Z burst model dis- 
cussed below is not an explanation of the highest energy cosmic rays. Detec- 
tion of GZK neutrinos could be essential for the understanding of very high 
energy cosmic ray protons. The normalization of the GZK flux is somewhat 
uncertain. Our results are based on a conservative value, that considers 
"standard" cosmological evolution For "strong" source evolution the 

same general behavior of the flux is valid, but the flux is a factor of 4 higher, 
leading to higher rates. Other studies [12] obtain fluxes different by up to 
one order of magnitude. Since the normalization we use is the smallest, our 
results are conservative and rates could actually be higher. 

In order to discuss how results depend on the shape of the initial flux, we 
also show two generic distributions, 1/E and l/E"^. The 1/E flux describes 
well the neutrinos from Z burst models ^3 El- Z burst models could 
explain the highest energy cosmic rays by extremely high energy neutrinos 



2 




Figure 1: Initial Neutrino Fluxes. 



scattering on nearby massive relic neutrinos. We also show results for the 
neutrino fluxes from AGN's for the model in Ref. 16J. This flux has an 
approximate behavior for E = 10^ — 10^ GeV, so towards the lower end 
of that energy range, our l/E"^ results apply. All of the fluxes used in our 
analysis are shown in Fig. Q except for the l/E"^ distribution, which would 
correspond to just a horizontal line. 

In the following section we describe our methods of computing the fluxes 
of neutrinos and charged leptons after propagation through the Earth. In 
section El we discuss the effect of regeneration and lepton energy loss on the 
propagation of the initial neutrino fluxes. In section |3] we present our results 
for the distributions of electromagnetic and hadronic showers from tau de- 
cays and neutrino interactions relevant to the different types of experiments. 
Other recent papers jSllIZI have focused on IceCube detection or on air show- 
ers. Related work has been recently presented at ICRC 2003 [TH]. Here, we 
consider IceCube, as well as RICE and ANITA, which also use the ice as 
detector and have higher energy sensitivities. In the case of ANITA, one has 
dramatically longer trajectories through the ice to consider. 
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2 Neutrino and Charged Lepton Propagation 



Neutrino attenuation and regeneration are governed by interaction lengths 
and decay lengths. In Fig. we show the neutrino interaction length (in 
water equivalent distance), as well as the tau decay length and the effective 
decay length when one includes tau energy loss in water jlHl EH] ■ 

The upper curve for the neutrino interaction length is equally applicable 
to antineutrinos, since at high energies the neutrino and antineutrino cross 
sections with nucleons are essentially equal because the cross sections are sea 
quark dominated. The neutral current contribution to the total cross section 
is about 1/2 of the charged current contribution. 

To compare the interaction lengths with physical distances, we note that 
D = 2R^cos0p ~ 6 X 10* cm.w.e. for 9 = 80° and D ~ 2 x 10'' cm.w.e. 
for 6 = 89° where = 6.37 x 10* cm is the radius of the Earth. Neutrino 
attenuation is clearly an important effect, even for nearly horizontal incident 
neutrinos. In the figures below, we mainly show results for a nadir angle 
of 85° where attenuation and regeneration effects in the propagation of u^-^s 
is in effect without dramatically reducing the flux. We also compare fluxes 
with those from incident 80° and 89° nadir angles. 

The effective decay length of the tau shows that, for energies above about 
10* GeV, the tau is more likely to interact electromagnetically than to decay 
|20j . For an initial tau energy of 10^^ GeV, the average energy just before it 
decays is a few x 10* GeV, depending on the density of the material the tau 
is passing through. Its effective decay length is of order 50 km in water. We 
use a density of p = 0.9 g/cm^ for ice. 

For tau neutrinos, we take into account the attenuation by charged cur- 
rent interactions, the shift in energy due to neutral current interactions and 
the regeneration from tau decay. For tau leptons we consider their produc- 
tion in charged current z/^ interactions, their decay, as well as electromagnetic 
energy loss. The tau neutrino and tau fluxes satisfy the following transport 
equations: 
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Figure 2: Neutrino interaction length (solid line) and tau effective decay 
length neglecting energy loss (dashed line) and including electromagnetic 
energy loss in water (solid line). 



-^ = a + ,3E. (3) 

Here F^^{E,X) = dN^JdE and Fr{E,X) = dN^/dE are the differential 
energy spectra of tau neutrinos and taus respectively, for lepton energy E, 
at a column depth X in the medium defined by 

X= p{L')dL'. (4) 
Jo 

We use the average density p{L) of the medium along the neutrino path, 
as given by the Preliminary Earth Model [21]. We note that Antarctica is 
covered by a sheet of ice with a few km thickness, so that some of the neutrino 
trajectories can go mostly or even entirely through ice rather than rock. For 
the neutrino interaction cross-sections we use CTEQ6 parton distribution 
functions [221 their extrapolation in the regions of interest to us P^ . 
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The decay length of the tau is \f^'^{E,X,9) = •ycTT-p. We use the decay 
distributions dn/dE written exphcitly in Ref. 4j, and we approximate the 
energy distribution of the neutrino interaction process with da/dy{E,y) ~ 
a{E)6{y — (y)), with y = [E — E')/E for incident neutrino energy E and 
outgoing lepton energy E' and {y) = 0.2. 

Equation Q describes the approximate energy loss of the tau. The 
parameter a is determined by the ionization energy loss and is negligible 
at high energy. The parameter (3 is due to radiative energy loss through 
bremsstrahlung, pair production and photonuclear scattering. The photonu- 
clear scattering becomes dominant at very high energies. As a first ap- 
proximation we use f3 = 0.8 x 10^^ cm^/g for the r energy loss. There 
is a negligible change (< a few %) if one uses a more realistic P{E) = 
(0.16 + 0.069 logio(^r/GeV)) x 10"^ cmVg ^ except at very high ener- 
gies {Er > 10^^ GeV), where the flux is extremely low. We use an energy 
loss parameter (3 = 7 x 10~^ cm^/g for the muons. 

It should be added here that for z/g the additional scattering on electrons 
should be considered because of the W boson resonance at 6.3 x 10^ GeV. 
However, this resonance is extremely narrow (narrower than the energy res- 
olution of the experiments) and it contributes a negligible amount to total 
rates 01111. 

In 121] it has been pointed out that the secondary fluxes of Ue and z/^ 
from the tau decays could also have a significant contribution for a flux of 
mono-energetic neutrinos. In it has been shown however that secondary 
neutrinos have negligible contribution for large nadir angles (i.e. for 6 > 60°) 
for generic initial neutrino fluxes, 1/E and l/E"^. The fluxes of secondary 
neutrinos are described by a transport equation similar to Eq. (Q), with the 
decay distribution dn/dE characteristic to each of the secondary neutrinos 
and we include them here, even though they are very small. 

3 Neutrino and Tau Fluxes 

lYiom. Eqs. (jll2l3|) for r's, and from suitably modified equations for muons, 
we evaluate the charged lepton fluxes at the end of the trajectory of a neutrino 
incident with nadir angle 9. 

For most astrophysical sources, the neutrinos are produced in pion decays, 
which determine the flavor ratio z/g : z/^ : z/,- to be 1 : 2 : 0. After propagation 
over very long distances, neutrino oscillations change this ratio to 1 : 1 : 1 
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because of the maximal ^ mixing. For the GZK flux, z/g and z/^ 
incident fluxes are different because of the additional contributions from v^. 
from neutron decay and from [i^ decays [12j. Because of this, the flavor 
ratio at Earth is affected by the full three flavor mixing and is different from 
1:1:1. Given fluxes at the source F° , F° and F° , the fluxes at Earth 
become: 

= Fl-\^n'2e,,{2Fl-Fl-Fl) (5) 
Fu, = F^r = \{Fl + Fl) + \,in'2e,,{2Fl~ F^- F'^^) (6) 

where 9i2 is the mixing angle relevant for solar neutrino oscillations. We 
have assumed that 6'23, the mixing angle relevant for atmospheric neutrino 
oscillations, is maximal and ^13 is very small, as shown by reactor experi- 
ments, as well as atmospheric and solar data. We use the initial GZK flux 
evaluated by Engel, Seckel and Stanev for standard evolution and we 
get the fluxes at Earth from the equations ^ and (jB)). The fluxes of and 
Vt- are still equal, due to the maximal ^ mixing. The main effect of 
the three flavor oscillations is to transform some of the low energy z/g's into 

and Vt. This could be useful since IceCube has very good sensitivity for 
detecting tracks, which are enhanced in this case. 

To illustrate the effect of the charged lepton lifetime and electromagnetic 
energy loss in propagation, we begin with muon and tau distributions. Fig. 
E] shows the tau and muon distributions for the GZK neutrinos at the end of 
the trajectory through the Earth with a nadir angle of 85°, with and without 
taking into account the energy loss of the charged lepton. As noted above, 
the incident z/^ and fluxes are equal. It is clear from here that the energy 
loss is extremely important and strongly limits the fluxes at very high energy. 
Including energy loss, the difference between the r and /i flux above 10^ GeV 
comes from the Vr pileup from interaction and regeneration, as well as to the 
difference in decay length and energy loss for taus and muons. Most of the 
charged leptons are produced in the last step of neutrino propagation. The 
muon flux below 10^ GeV dominates the tau flux because of the much longer 
decay length of the muon. 

Fig. m shows the distribution of z/^, z/^ and z/g fluxes after propagation 
through the Earth for the same initial flux and angle. For z/g and z/^ the 
transport equations are effectively decoupled from those of the charged lep- 
tons and their propagation is given by Eq. (Q) in which the last term is not 
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Figure 3: Charged lepton distributions for GZK neutrinos, at a nadir angle 
of 85°. 

present. On the same figure we show the flux of secondary z/^ neutrinos, 
which is the same with that of the secondary v^.. It can be seen that their 
contribution is neghgible, with maximum corrections of the order of a few 
percent around 10^ GeV and much smaller for most energies. At larger nadir 
angles, the secondary neutrino contribution will be an even smaller fraction 
of the z/g and fluxes. At smaller nadir angles, the secondary flux will be a 
larger percentage of the overall neutrino flux, but the flux will become more 
strongly attenuated. 

Fig. El shows the same distribution as in the previous figure, together 
with the distributions obtained if energy loss of tau leptons or the shift in 
energy due to neutral current is not included in the propagation. Electro- 
magnetic energy loss is extremely important at energies above 10^ GeV and 
this effect causes strong suppression of the neutrino fluxes at very high en- 
ergy. If the tau would not lose energy in its propagation, the regeneration 
effect would be much bigger at high energy. The i/^ — > r — » z/^ regeneration 
is a significant effect for energies between 10^ and 10^ GeV, as can be seen 
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Figure 4: Neutrino distributions for GZK neutrinos, at a nadir angle of 85°. 



more clearly in Fig. [7| The neutral current interactions of the neutrinos are 
also very important, as can be seen in Fig. El by comparing the flux of z/^ af- 
ter correct propagation with a simple attenuation with exp{-~D / C\^^[E)) for 
column depth D and interaction length C evaluated using the charged cur- 
rent cross section. The difference can be as large as two orders of magnitude 
at very high energy. 

Fig. ini shows the r and Vr distributions for the same initial neutrino flux, 
but for different nadir angles. As the neutrinos pass through more matter, 
the regeneration effects become relatively more important, but the overall 
fluxes get significantly attenuated. At a nadir angle of 80° the Earth is 
already opaque to neutrinos with energies above 10^° GeV. 

Fig. [7| shows the ratio of the tau neutrino flux after propagation to 
incident tau neutrino flux, for 89°, 85°and 80°. This ratio illustrates a com- 
bination of the regeneration of v^. due to tau decay and the attenuation of all 
neutrino fluxes. For 89°, where both the total distance and the density are 
smaller, the attenuation is less dramatic, and the flux can be significant even 
at high energy. The regeneration in this case can add about 25% corrections 
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Figure 5: Tau neutrino distributions for the GZK flux, at a nadir angle of 
85°. 



at energies between 10^ and 10^ GeV. For 85° the relative enhancement is 
around 80% and peaked at slightly lower energies, while at 80° it is almost a 
factor of 3 at low energy. At 80°, however, the flux is very strongly attenuated 
for energies above a few xlO'' GeV. 

It is already clear from Fig. [7|that the total rates will be dominated by 
the nearly horizontal trajectories that go through a small amount of matter. 
The largest pileups occur when the trajectory of the neutrino passes through 
a larger column depth. For the higher energies relevant to RICE and ANITA, 
this doesn't translate to higher fluxes of z/t-'s. Attenuation is the main effect 
at those energies. Rates can get significant enhancements at low energies 
where the regeneration from tau decays adds an important contribution even 
for longer trajectories. 
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Figure 6: The r and Vr distributions for GZK neutrinos for different nadir 
angles. 

4 Showers 

The above distributions illustrate the effects of the propagation through mat- 
ter on the initial neutrino fluxes. We now turn to the experimentally relevant 
observables, which are the showers. 

We consider hadronic showers produced in interactions, as well as 
electromagnetic and hadronic showers from r decays in ice. All other flavors 
of neutrinos will contribute to the hadronic showers, while Ve will also con- 
tribute to electromagnetic showers. The shower distributions will be different 
for each type of experiment, depending on the detector setup. 

We start by considering contributions to electromagnetic and hadronic 
showers in a depth of ~ 1 km of ice, relevant for experiments like IceCube 
and RICE. 

We use Eqs. (jUEl) and the corresponding ones for e and flavors to obtain 
the tau and neutrino distributions {F^-^Et-) and F^{E^)) after propagation up 
to the last few km of ice. We then estimate the shower distributions that 
are produced over the distance d ~ few km. The contributions of neutral 
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Figure 7: Ratio i/finay ^initial f^^. qiK neutrinos, at nadir angles of 80°, 85° 
and 89°. 

currents and regeneration over a few kilometers is small and here we neglect 
them for the last part of the trajectory, which is where the showering takes 
place. This can no longer be done when showers are detected over long 
distances, as will be the case for ANITA. 

There are two contributions to electromagnetic showers from tau decays. 
One is from the decay of the taus produced outside the detector that decay 
electromagnetically inside the detector, giving "lollipop" events in IceCube 
and electromagnetic showers in RICE. The distribution of these showers is 
given as a function of the energy of the electron produced in the decays by: 



where F{E) = dN/dE is the flux of taus just before the final column depth 
d in ice and p = 0.9 g/cm^ is the density of ice. For negligible energy loss 
(which is a very good approximation below 10^ GeV in 1 km of ice) this 
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becomes: 

FZ'm = dE^FAEr) (l - e-'^^/^^-^^^)) ^(E., E,) (8) 

The other contribution comes from the taus that are produced and decay elec- 
tromagnetically inside the detector. In IceCube these correspond to "double 
bang" events. The distribution is obtained from: 

rdp j-dp l-Emax l-Ev 

FX;\E,) = dx dy dE,l dE,F,^{E,^) 

Jo Jx JEe JEe 



1 (in^_ 



Xf^{E,e-Py) dE^ 
When energy loss is neghgible, this becomes: 



[E.^ErC-fy) (9) 



FZ'iE^) = j^^ dE^J^^ dE^F^AE^^)N^^{E^^,E^)^{E^,E^ 



The main source of electromagnetic showers comes from the z/g charged 
current interactions given by: 



rEmax , nr,^ ^ 1 d.rr^'-' 



Hadronic showers from tau decays are given by Eqs.(jHI) with the electron 
replaced by hadrons. 

The hadronic showers from i/j. interactions are given by: 

F^H,. = dEMK){l - e-''"'^''^''^^)^,^{E., En). (12) 

J Eh a{hu)<ii^h 

In Eq. fll2p . a is the total cross section and the hadronic energy comes from 
the energy transfer to the target nucleus. 
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Hadronic showers are also produced in z/^ and neutral current and 
charged current interactions. These are also obtained from Eq. (jl2p . using 
the corresponding fluxes. 

We briefly describe below all types of interactions and their signatures in 
the various experiments. We then show shower distributions relevant in each 
case. 

4.1 Showers in kilometer size detectors: IceCube, RICE 

A km'^ Cherenkov detector in ice, IceCube [7j has sensitivity to high energy 
neutrinos up to ~ 10* — 10^ GeV, the energy range where distinguishing the 
effects of tau neutrinos from muon and electron neutrinos could be possible. 

Neutral current interactions of all types of neutrinos produce hadronic 
showers. Charged current interactions of Ve contribute to both hadronic 
and electromagnetic showers. It is likely that IceCube would observe a big 
shower with the total energy of the incoming neutrino rather than be able 
to separate the two components (hadronic and electromagnetic). Charged 
current interactions of give a hadronic shower and a muon. The muon 
track is the main signal in IceCube, so these events are easily identifiable, 
containing a shower and a long track emerging from it. At high energies 
however, the same signal could be produced by a Vj. interaction, since the 
decay length of the tau becomes longer than the size of the detector for 
energies above a fewxlO^ GeV. 

Charged current Vr interactions can have very different signatures depend- 
ing on energy. At 10^ GeV, the tau decay length is ~50 m, and the shower 
(hadronic or electromagnetic) from the tau decay cannot be separated from 
the hadronic shower from the initial interaction. At a few times 10^ GeV 
the range of the tau becomes a few hundred meters and can give the charac- 
teristic signal of the "double bang" events |23]. These are events where the 
shower from the neutrino interaction and the shower from the tau decay can 
be separated and are both observed in the detector, together with the tau 
track. At a fewxlO^ GeV the tau decay length is already longer than 1 km, 
and the taus look like muons. What is seen is the shower from the neutrino 
interaction and then a track. 

Taus produced in neutrino interactions outside the detector can generate 
"lollipop" events. The initial shower from the neutrino interaction that pro- 
duces the tau is missed and what is seen is the track of the tau and the shower 
from its decay that IceCube can identify as electromagnetic or hadronic from 
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the existence of muon tracks in the shower. 

RICE jH] uses dipole antennas in the Antarctic ice to measure radio fre- 
quency Cherenkov radiation from high energy showers. For energies around 
10^ GeV the effective volume of the detector is ~ 15 km^sr. At present, 
RICE limits the fluxes of downward Ve with energies between few x 10^ and 
10^^ GeV jH]. Interactions modeled by a Monte Carlo at an energy of 10^ GeV 
are detected from within a depth of 1 km from the surface and out to a radial 
distance of about 4 km. The effective volume of the detector has a strong 
dependence on energy below 10^ GeV. Hadronic and electromagnetic show- 
ers can be separately identified, with somewhat different effective volumes, 
depending on energy. The RICE experiment can therefore measure electro- 
magnetic showers from v^, interactions, hadronic showers from the charged 
and neutral current interactions of all flavors of neutrinos and electromag- 
netic and hadronic showers from tau decays. 

Fig. IHl shows the electromagnetic shower distributions at a nadir angle of 
85°. In the absence of the tau neutrinos, only z/g interactions lead to elec- 
tromagnetic showers. These i^e ^ e CC conversions still dominate at high 
energies, even when tau neutrinos produce tau leptons that decay electromag- 
netically. At lower energies however, the taus add a significant contribution 
to the electromagnetic shower signal. Tau decays give their most impor- 
tant relative contribution to electromagnetic showers at electron energies of 
a fewxlO''GeV. These decays are separated in "double bang" and "lollipop" 
events on the same plot. It can be seen that the contribution of the "double 
bang" is relatively small. They do not contribute significantly to total event 
rates, but even in very small numbers they are important as a characteristic 
signature for taus. Secondary z/g interactions also produce electromagnetic 
showers, but it can be seen that their contribution is very small. 

Fig. ini shows the hadronic shower distributions. The features seen in the 
neutrino fluxes are recovered in the hadronic showers as well. The shower 
rates from z/,-, z/^ and Ve are nearly identical at shower energies above 10^ 
GeV. The tau neutrino pileup is important below 10^ GeV. Between 10^ 
and 10^ GeV tau decays give the main contribution to the hadronic shower 
rate. Showers from z/e and z/^ secondary neutrinos are also shown and their 
contribution is negligible. 

Fig. ITUl shows the ratio of the electromagnetic shower rates at nadir angle 
85° in the presence and absence of oscillations for the GZK and Z burst neu- 
trino fluxes. This ratio illustrates the effect of the oscillations on the signal 
and in particular the possible enhancements due to tau pileup. In absence 
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Figure 8: Electromagnetic shower distributions for GZK neutrinos, at a nadir 
angle of 85° for a km size detector. The r decay curve is the sum of the 
lollipop and double bang curves. 



of oscillations, the only contribution to electromagnetic showers comes from 
i/g interactions. In the presence of — > Ur oscillations, electromagnetic de- 
cays of taus from tau neutrinos add significant contributions to these rates 
at energies below 10*^ GeV. In the same time, for the GZK flux, Ve z/^,r 
oscillations reduce the number of z/g's at low energy, such that below a few 
xlO^ GeV there are fewer electromagnetic showers than in the absence of 
oscillations. 

Fig. im shows the ratio of the hadronic shower distributions in presence 
or absence of oscillations for the same fluxes and nadir angle. Without oscil- 
lations, the hadronic showers come from z/g and interactions. In presence 
of oscillations and the z/g and secondaries from tau decay contribute as 
described above. For hadronic showers the enhancement due to the presence 
of Vr is ~ 50% at energies of a few xlO^ GeV. For hadrons there is also some 
small enhancement for shower energies between 10^ and 10^ GeV, which is 
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Figure 9: Hadronic shower distributions for GZK neutrinos, at a nadir angle 
of 85° for a km size detector. 



not present for the electromagnetic showers. 

The tau contribution can be even more clearly seen if some of the other 
flavors can be separately identified. As previously discussed, IceCube can 
identify showers from i/^ and z/^ charged current interactions by the /i or r 
tracks that emerge and exit the detector. Once these can be removed, what 
remains are hadronic showers from neutral current interactions of all flavors 
of neutrinos, v^. CC interactions in the detector, as well as the hadronic 
showers from r decays in the detector (both from taus produced within the 
detector and from taus produced outside the detector that propagate in and 
then decay). The enhancement in the ratio is even higher in this case, about 
80%. 

Fig. ^1 also shows the hadronic shower distributions for GZK neutrinos, 
but at a nadir angle of 89°. As expected, there are more showers for this 
distance than for the longer one, but the relative contribution of the tau is 
smaller in this case. This is because both the attenuation and the regenera- 
tion effects are much smaller for smaller column densities. 
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Figure 10: Ratio of electromagnetic shower rates in the presence and absence 
of z/^ Uj. oscillations for GZK and 1/E neutrino spectra for nadir angle 
85° for a km size detector. 



The experimental angular resolution for high energy showers is about 10° 
to 25°. Averaging over such angles the effect of the pileup is somewhat re- 
duced (to ~ 10% of the event rate), since total rates are dominated by the 
trajectories that go through a small amount of matter. However, rates can 
get significant enhancements at low energies where the regeneration from 
tau decays adds an important contribution even for longer trajectories. Con- 
sequently, an experiment with low energy threshold has a better chance of 
detecting the effect. IceCube has a low energy threshold and an energy res- 
olution of ~ 10%, so it is in a good position to look for pileup effects, given 
a high enough neutrino flux. 

The same general features remain true for different initial fluxes. How- 
ever, the energy distributions are different for each and the enhancement due 
to regeneration appears at somewhat different energies. The regeneration 
effects are also smaller for the fluxes steeper at high energy. For example, 
from Figs. ITUl and ITD it can be seen that for the 1/E flux, predicted, for 
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Figure 11: Ratio of hadronic shower rates in the presence and absence of 
Vj. oscillations for GZK and 1/E fluxes for nadir angle 85° for a km 
size detector. 



example, for Z burst (ZB) models [H], the pileup is slightly bigger, but at 
a lower energy than for the GZK neutrino flux. This is just as expected, 
because the GZK flux is steeper at higher energy, but does not fall as rapidly 
at intermediate energies. 

Fig. El shows the electromagnetic shower distributions for a.1/ E flux, for 
example predicted for Z burst and for a. 1/ E"^ generic flux for a nadir angle 
of 85°. The l/E"^ distribution has a normalization of 10^*^ GeV~^ km~^ yr~^ 
sr~^, about an order of magnitude below the present AMANDA limit j^ . 
For the ZB flux the normalization is 1 GeV~^ km~^ yr~^ sr up to 2.5 x 10^^ 
GeV. Above this energy the flux is cut off and drops as 1/E^. The steepness 
of the l/E"^ flux results in a small pileup and thus the relative contribution of 
taus is much smaller in this case. Fig. IT^ shows the same shower distributions 
for a nadir angle of 89°. Like for the GZK flux, the attenuation in this case 
is smaller than for 85°, particularly at high energy. The effect of attenuation 
at high energy is striking especially for the Z burst flux. At 85° there is 
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Figure 12: Hadronic shower distributions for GZK neutrinos, at a nadir angle 
of 89° for a km size detector. 

almost no flux left at energies above 10^° GeV, while at 89° the neutrino 
flux is almost unattenuated and still has the 1/E shape, being orders of 
magnitude higher than for the longer pathlength. However, the effects of 
regeneration become smaller for 89° even for the less steep Z burst flux. It 
can be clearly seen that, depending on the energy threshold of the detectors, 
the contribution to event rates comes from different trajectories. At high 
energies the paths that go through the material with small column density 
will dominate event rates. However, at energies below 10^ GeV trajectories 
that go through 10-15 times more material can contribute equally due to the 
pileup and additional tau decays. 

Fig. El shows the hadronic shower distributions for the fluxes used in 
the previous figures, for a nadir angle of 85°. The energy distribution of the 
hadronic showers is different from that of electromagnetic ones due to the dif- 
ferent decay distributions. Taus decay mostly to hadrons and these hadrons 
carry most of the energy of the taus. Consequently, the tau contribution to 
hadronic showers is much higher, bigger than the z/^ contribution at energies 
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Figure 13: Electromagnetic shower distributions for nadir angle 85° for 1/ E 
and 1/ E"^ characteristic fluxes for a km size detector from Ve interactions and 
from r decays. 



around a few x 10^ GeV, and extends to higher energies. 

Fig. ^1 shows the electromagnetic distributions for the AGN flux in |Tn] 
for nadir angles of 85° and 89°. This flux is also steep at high energy and 
consequently the regeneration effects are small. The relative contribution of 
the taus is bigger for 85°, but the overall rates are higher for 89°. 

The energy threshold of RICE is high and for the energies where the 
experiment has good sensitivity the Vr enhancemnt is limited by the long 
lifetime of the tau lepton and by its energy loss, such that it cannot be 
observed. We want to investigate if this remains true in the case of ANITA, 
which has a much larger effective volume and consequently has the potential 
to detect many more tau's decaying over long trajectories. 
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Figure 14: Electromagnetic shower distributions for nadir angle 89° for 1/ E 
and l/E"^ characteristic fluxes for a km size detector. 

4.2 Showers in ANITA 

The ANITA experiment also uses the ICG cLS cL neutrino converter ^U]. The 
long duration balloon missions will monitor the ice sheet from 40 km in 
altitude to a horizon approaching 700 km for refracted radio frequency signals 
with an effective telescope area of IM km^. The geometry of the experiment 
is rather complicated, as it has to take into account the Cherenkov angle 
of the radio emission with respect to the particle trajectory, its refraction 
at the ice surface and the position of the balloon. The huge volume covered 
gives ANITA remarkable sensitivity for detecting very high energy neutrinos. 
While IceCube and RICE can only detect showers produced in a km^ volume, 
ANITA can detect all the showers produced over long distances in ice. 

For example, for a nadir angle of 89°, the entire trajectory (~ 222 km) of 
the neutrino is in ice, at less than 1 km depth, and all showers produced over 
this distance can be observed. For a nadir angle of 85°, observable showers 
(at less than 1 km depth) could be produced over a distance of ~12 km. 



22 



Hadronic showers 



le+08 



T, 1/E 



[e+06 le+07 le+08 le+09 le+10 le+11 

E [GeV] 




le+12 



Figure 15: Hadronic shower distributions for nadir angle 85° for 1/E and 
l/E"^ characteristic fluxes for a km size detector from neutrino interactions 
and tau decays. The z/g distributions are identical to the distributions. 
For the 1/ E'^ flux all neutrinos distributions are almost the same. 



Given the very large detector trajectory is not fully defined by the 

nadir angle of the incident neutrino at the entrance point, one also needs 
the position of the entrance/exit point with respect to the icecap and the 
balloon. The trajectories that maximize the path through ice, rather than 
the ones that combine ice and rock, are the ones likely to dominate event rates 
because of the much smaller attenuation. Those going through more rock give 
larger relative enhancements due to regeneration. They could in principle 
contribute as much as the others in the region where the regeneration is 
effective. 

As previously discussed, due to the showering over entire trajectories, we 
can no longer use Eq. (jSHT^ , but rather we have to combine the propagation 
and showering from the beginning in order to correctly take into account neu- 
tral current interactions, energy loss of tau leptons and neutrino regeneration 
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Figure 16: Electromagnetic shower distributions for the MPR AGN model 
[T^ for a km size detector. 

from tau decay. 

Fig. El shows the electromagnetic shower distributions for the GZK and 
ZB initial fluxes over a trajectory of 222 km in ice. Qualitatively, these 
are similar to the showers obtained in kilometer size detectors for the same 
trajectory. However, in this case the rates are larger by up to three orders 
of magnitude, depending on energy, due to the much larger surface of the 
detector. The energy distribution is somewhat different, being more spread 
out toward higher energies. Due to the longer detectable pathlength, higher 
energy taus, with decay length much longer than 1 km, can now decay pro- 
ducing observable showers. 

Fig. ITHl shows the hadronic shower distributions for the same initial fluxes 
and trajectory. The regeneration of the neutrinos in this case is very small 
and distributions are almost the same as the z/^ ones. For the ZB flux, 
the Ve distribution is the same as for z/^, while for the GZK neutrinos it is 
different due to the difference in the initial fluxes. 
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Figure 17: Electromagnetic shower distributions for detection over 222 km 
of ice. 

Fig. El shows the ratio of electromagnetic and hadronic shower rates 
in the presence and absence of — >■ Vr oscillations corresponding to the 
distributions in the previous figures. The maximum enhancement due to the 
presence of Vr is about 40% at this angle, as expected since this trajectory 
has low column density. However, as previously discussed, the enhancement 
occurs in a larger energy range and it peaks at higher energy than in the case 
of a small size detectors. In 1 km only taus with energies below a few xlO^ 
GeV have a significant probability to decay, while much higher energy taus 
can decay over the total distance of more than 200 km. 

5 Conclusions 

We have studied in detail the propagation of all flavors of neutrinos with 
very high energy {E > 10^ GeV) as they traverse the Earth. Because of the 
high energies, we have limited our consideration to nadir angles larger than 
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Figure 18: Hadronic shower distributions for detection over 222 km of ice. 
The distribution for GZK neutrinos is the same with that of Ur. For the 
1/E flux all neutrino distributions are almost the same. 

80°. We are particularly interested in the contribution from tau neutrinos, 
produced in oscillations of extragalactic muon neutrinos as they travel large 
astrophysical distances. After propagation over very long distances, neutrino 
oscillations change an initial (source) flavor ratio of 1 : 2 : to 1 : 1 : 1 
because of the maximal z/^ mixing. For GZK neutrinos, the flavor ratio 
at Earth deviates from 1:1:1 because the incident fluxes are different. At lower 
energies and smaller nadir angles, tau neutrino pileups from regeneration 
via — > T ^ z^r E] euhauced electromagnetic and hadronic signals in 
kilometer-sized detectors. Our aim here was to see if there are similar effects 
at high energy. 

In our propagation of neutrinos and charged leptons through the Earth, 
we have focused on kilometer-sized neutrino detectors, such as ICECUBE 
and RICE and on a detector with much larger effective area which uses 
Antarctic ice as a converter, ANITA. Our study can easily be generalized to 
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Figure 19: Ratio of electromagnetic and hadronic shower rates in the presence 
and absence of z/^ oscillations for GZK and 1/E fluxes of neutrinos for 
detection over 222 km of ice. 

other experiments and to propagation in materials other than ice. 

We have found that the flux above 10^ GeV resembles the flux. 
The lore that the Earth is transparent to tau neutrinos is not applicable in 
the high energy regime. Tau neutrino pileups at small angles with respect to 
the horizon are significantly damped due to tau electromagnetic energy loss 
above Er ~ 10^ GeV if the column depth is at least as large as the neutrino 
interaction length. 

At lower energies, E < 10^ GeV, regeneration of becomes important 
for trajectories where the other flavors of neutrinos are strongly attenuated. 
The regeneration effect depends strongly on the shape of the initial flux and 
it is larger for flatter fluxes. The enhancement due to regeneration also 
depends on the amount of material traversed by neutrinos and leptons, i.e. 
on nadir angle. For GZK neutrinos, we have found that the enhancement 
peaks between 10^ and a fewxlO'' GeV depending on trajectory. For 85° the 
enhancement is about a factor of two, while for 80° it is a factor of three. 

We have translated the neutrino fluxes and tau lepton fluxes into rates 
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for electromagnetic and hadronic showers at selected angles to see the effect 
of attenuation, regeneration, and the different energy dependences of the 
incident fluxes. We have focused on comparing the v^- contribution to the Ve 
and contributions to determine in what range, if any, t/^'s enhance shower 
rates. 

The Vt flux enhancements depend on the shape of the initial flux. The 
electromagnetic showers are more sensitive to this shape than hadronic ones. 
The relative enhancement in hadronic showers is also smaller than for the 
electromagnetic showers. This is because for the electromagnetic signal the 
only contribution in the absence of taus is from electron neutrinos, while for 
hadrons the tau contribution is compared to a much larger signal, from the 
interactions of all flavors of neutrinos. We have included contribution from 
secondary neutrinos, which we find to be relatively small for all fluxes. 

For kilometer-sized detectors, at for example a nadir angle of 85°, the 
maximal enhancemnet due to v-^ contribution to electromagnetic shower rates 
for the GZK flux is about 50% at 3 x 10^ GeV, while for the l/E flux, it 
is even larger, about 70%, at slightly lower energy. In the case of hadronic 
showers for which the events identified by muon tracks have been removed, 
the contribution peaks at about 2 x 10^ GeV and it gives an enhancement 
of about a factor of 1.8 for the GZK flux. These energy ranges are relevant 
for IceCube, but not for RICE. For energies relevant to RICE, tau neutrinos 
do not offer any appreciable gain in electromagnetic shower signals compared 
to z/e ~^ e CC interactions, and they contribute at essentially the same level 
as to hadronic shower rates through NC interactions. 

One of the reasons that tau neutrinos do not contribute large signals to 
kilometer-sized detectors at very high energies is that high energy tau decay 
lengths are very large, so the probability of a tau decaying in the detector is 
low. For detectors like ANITA which can sample long trajectories through 
the ice one would expect a larger tau neutrino contribution to the signal from 
tau decay. Despite the long trajectory (222 km with a maximum depth of 1 
km for a neutrino incident at 89° nadir angle) the tau contributions to the 
electromagnetic shower rate is quite small for fluxes expected to contribute 
in the ANITA signal. For hadronic showers, the suppression of r decay to 
hadrons relative to NC interaction contributions is about the same as for 
electromagnetic showers compared to Ve e. The contribution to the 
hadronic shower rate from interactions is the same as the Ve contribution. In 
summary, for ANITA, tau neutrinos do not give any additional signal beyond 
what one would evaluate based on no regeneration from i/^ — > r — > i^,- due to 
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tau electromagnetic energy loss at E > 10^ GeV. 

In addition to the experiments discussed here, there are many studies 
concerning the possibility for detection of radio Cherenkov emission from 
showers in materials other than ice. It has been noted [27j that rock salt 
formations have similar properties to the Antarctic ice and can therefore be 
used as large scale neutrino detectors. Salt has a higher density {psait = 2.2 
g/cm^) than ice (pjce = 0.9 g/cm^), so it is possible to achieve an effective 
detection volume of several hundred km'^ water equivalent in salt. This is 
somewhat larger than RICE, achieved with a much smaller actual detector 
size. The threshold for detecting the radio signal from showers in salt is of 
the order of ~ 10^ GeV, similar to RICE, but lower, such that detection of 
extra signals from enhancements would be more promising. 

Also proposed is LOFAR [2H1, a digital telescope array designed to detect 
radio Cherenkov emission in air showers. LOFAR has sensitivity in an energy 
range of ~ 10^ — 10^^ GeV, so it can detect showers at much lower energies 
than other radio Cherenkov experiments. LOFAR will likely be configured 
to detect horizontal showers from skimming neutrinos as well. With its low 
energy threshold, LOFAR has an excellent opportunity to observe the shower 
enhancement at lower energies due to Ur regeneration and tau pileup, which 
is not easily accessible in ANITA. 
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